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Abstract 

We generalize a previously known class of exact an- 
alytic solutions of relativistic perfect fluid hydrody- 
namics [1] for the first time to arbitrary temperature- 
dependent Equation of State. We investigate spe- 
cial cases of this class of solutions, in particular, we 
present hydrodynamical solutions with the Equation 
of State determined from lattice QCD calculations. 
We discuss the phenomenological relevance of these 
solutions as well. 

1 Introduction 

The interest in relativistic hydrodynamics grew in 
past years mainly due to the discovery of the almost 
perfect fluidity of the experimentally created Quark- 
Gluon-Plasma [2]. Hydrodynamical models aim to 
describe the space-time picture of heavy- ion collisions 
and infer the relation between experimental observ- 
ables and the initial conditions. Besides numerical 
simulations there is also interest in models where ex- 
act solutions of the hydrodynamical equations are 
used. 

In this paper we generalize a previously known 
class of exact solutions of relativistic perfect fluid hy- 
drodynamics [1] to the case of arbitrary, temperature 
dependent speed of sound. The mentioned class of so- 
lutions form the basis of the relativistic Buda-Lund 
hydrodynamical model [3]. This model yields a suc- 
cessful description of hadronic observables at RHIC 
energies (such as the pseudorapidity and transverse 
momentum dependence of the azimuthal anisotropy 
of different hadrons as well as the HBT radii [3]), 
and the reconstructed final state in this model cor- 
responds to simple explicit scaling solutions of hy- 
drodynamics. The same final state however can be 
achieved from many initial states, depending on the 
Equation of State [4]. If one is given a temperature 



dependent speed of sound as Equation of State, the 
solution presented in this paper thus can be used to 
determine the initial state from the reconstructed fi- 
nal state of a heavy-ion collision. As an example, 
we describe the time dependence of the system if one 
assumes the Equation of State from lattice QCD. 

The solutions given in this paper are the first ex- 
act analytic solutions of 1+3 dimensional relativistic 
hydrodynamics, to utilize an arbitrary Equation of 
State. 1 



2 Basic equations 

Let us adopt the following notational conventions: 
the fluid coordinates are x M = (t, r), where r = 
{ r x, f y , r z ) is the spatial coordinate, and the met- 
ric tensor is — diag (1, — 1, — 1, — 1). (We de- 
note space-time indices by Greek letters, space in- 
dices by Latin letters, and assume the summation 
convention.) The fluid four-velocity is = 7(l,v), 
where v is the three- velocity and 7 = 1/y/l — v^. 
The thermodynamical quantities are denoted as fol- 
lows: p is the pressure, e is the energy density, a is 
the entropy density, T is the temperature. If the fluid 
consists of individual conserved particles, or if there 
is some conserved charge, then the conserved number 
density is denoted by n, and the corresponding chem- 
ical potential by fi. (For more than one conserved 
number densities, we may use indices to distinguish 
them.) All these quantities have dependence on x M , 
but mostly this will not be explicitly written out. 

The basic hydrodynamical equations are the conti- 



1 Note that it has been discussed in Ref. [5] that the en- 
tropy flow can be calculated with an arbitrary EoS (speed of 
sound) from the Khalatnikov-potential, once the solution of 
the general Khalatnikov equation is known. 
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nuity and energy-momentum-conservation equations: 

(nu") = 0, (1) 
d v T^ v = 0. (2) 

The energy-momentum tensor of a perfect fluid is 

= (e+p) u^u v - pg^. (3) 

Eq. (2) can be then transformed to (by projecting it 
orthogonal and parallel to u^, respectively): 

(e + p) u v d v u» = W - uV) d vP , (4) 
{e+p)d v u v + u v d v e = Q. (5) 

Eq. (4) is the relativistic Euler equation, while Eq. 
(5) is the relativistic form of the energy conservation 
equation. In Appendix A we recall the well-known 
fact that Eq. (5) is equivalent to the entropy conser- 
vation equation: 



d p (tro") 



0. 



(6) 



An analytic hydrodynamical solution is a func- 
tional form of e, p, T, (and, if dealt with, n), 
which solves Eqs. (4) and (5), and, if present, n also 
solves Eq. (1). The quantities e, p, T, and also a, 
and n are subject to the Equation of State (EoS), 
which closes the set of equations. We investigate the 
following EoS: 



e = n(T)p, 



(7) 



and for the case when there is a conserved n number 
density, the additional assumption is 



p = nT. 



(8) 



For the case of k (T) = n constant, an ellipsoidally 
symmetric solution of the hydrodynamical equations 
is presented in Ref. [1]: 



with Xq, Yq, Zq constants. The quantity s has el- 
lipsoidal level surfaces, and obeys u v d v s = 0. We 
call s a scaling variable, and V the effective volume 
of a characteristic ellipsoid 2 . This solution is non- 
accelerating, ie. obeys u^d^u^ = 0. In the next sec- 
tion we present a generalization of this class of solu- 
tions to more general EoS. The new solutions will be 
presented in Section 4, while Section 3 details their 
derivation. 



3 General Equation of State 

In order to find more general solutions, where a tem- 
perature dependent EoS can be used (as in Eq. (7)), 
for a given w M velocity field we may define the V and 
s quantities by their properties that 



0. 



(14) 



With these quantities, Eq. (1) is automatically solved 
(for the case when there is a conserved charge 
present) if 



V 

n —v{s) 



(15) 



again, with arbitrary v (s) function. To solve the (5) 
energy equation, we must make a distinction between 
two possible cases. The first case is if we take a con- 
served n into account, and use the EoS e = k (T)p, 
p = nT as in Eqs. (7) and (8). The second case is 
when we do not consider any conserved n. In Ap- 
pendix A we show that in both of these two cases 
the energy equation Eq. (5) can be transformed to 
an equation for T: in the first case with conserved n, 
we have 



Al = ^, t = Vt 2 - r 2 = (9) 



d(«T) d^T 



AT 



V 



0. 



(16) 



n = n a ^v{s), T = T (^j (10) 
with v (s) being an arbitrary function and 

S=4^ + T^ + ^, V = T\ (11) 



X 2 Y 2 



Z 2 



where X, Y, and Z are the time (t) dependent prin- 
cipal axes of an expanding ellipsoid. They have the 
explicit time dependence as 

X = X t, Y = Y t, Z = Z t (12) 



while in the case where there is no conserved n, we 
have 



V 



1 dn 



K 



1 dT T 



= 0. 



(17) 



2 Note that with the X (t), Y (t), Z (t) time-dependent axes 
introduced as here, we can write the velocity field as 



Y 
— i 

Y 



(13) 



which underlines the resemblance of this solution to certain 
non-relativistic exact solutions with Hubblc-likc expansion [6] . 
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Remarkably, these equations are not the same (how- 
ever, we may note that in the case when k = 
const., they yield the same condition). We call these 
cqutions the temperature equations for the two cases. 
With the introduction of the / (T) function as 



f(T) = exp 




l_d_ 

£d£ 



M/3) P])dp 



(18) 



for the case of conserved n, and as 



f(T) = exp 




d«(/3) 



k(P) + 1 d/3 



d/3 



(19) 



for the case of vanishing n, the temperature equations 
can be cast in the following universal form: 



/CO 



V 



= 0. 



(20) 



For any given k(T) function we can determine f{T), 
and write up the solution of the above equation as 



t = r 1 



Us) 



(21) 



with arbitary £ (s) function. (For convenience, we 
may normalize £ (s) so that £(0) = 1.) Knowing that 
u^d^s — 0, it is easy to see that this indeed solves 
Eq. (20). Note that if k = const., then due to Eq. 
(21): 



/(T,= [ % 



(22) 



As a generalization of the solution recalled in the 
previous section, we assume that and thus s and 
V has the same forms as in Eqs. (9) and (11). 



and using the expression of T from Eq. (21) and the 
definition of and V from Eq. (23), we find that it 
is equivalent (for any k (T), thus for any / (T)) to 



(26) 



In the case of non- vanishing n, using Eq. (15) and 
p = nT, the Euler equation for our non-accelerating 
velocity field transforms to the following equation: 

Td^n + nd^T = Tu^u v d v n + nu^u u d v T. (27) 

Substituting n and T from Eqs. (15) and (21), and 
the definition of V, we get from this equation the 
following constraint: 



d^s = 0, 



(28) 



where we have introduced the following function: 

yf^'iy) 



(29) 



Since d^s ^ 0, we see from Eq. (28) that there are 
two cases: for any EoS (ie. for any k (T) and thus any 
ip(T) function) we get a solution if v (s) = £ (s) = 1. 
The other possibility is if k = const. It is easy to see 
that this case is equivalent to ip = re -1 = const, and 
so Eq. (28) is solved if £ = V - X I K and so from Eq. (22) 
we get T = T (Vq/V) 1 ^ ^ _1 (s), i.e. the same as in 
Eq. (10). In this case we indeed obtain the known 
solution of Ref. [1], recited in Eqs. (9)-(ll). 



— , V 



Xgt* 
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Y 2 t 2 



+ 



(23) 



For this velocity field, u^d^u^ = 0, so the remaining 
equation, the Euler equation of (4) is equivalent to 



d^p = u fl u"d u p. 



(24) 



In the case of vanishing n, using the thermodynamic 
relation dp — adT, Eq. (24) simplifies to 



(25) 



New solutions for 
Equation of State 



general 



Summarizing and rewriting the results presented in 
the previous section, we found new solutions to the 
relativistic hydrodynamical equations for arbitary 
e = k (T)p Equation of State, and these are the first 
solutions of their kind (i.e. with a non-constant EoS). 
In the case where we do not consider any conserved 
n density, the solution can be presented in the fol- 
lowing form, in terms of a and T, with T given 
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in an implicit form: 



T 



exp 




k 08) + 1 d/3 



d/3 



b 5 
y 



exp 




CTo- 



— =exp 



K (/3) , d«08) 







d/3 



d/3 y K 08) + 1 



Also, for the case when the pressure is expressed as 
p = nT with some conserved n density, the new so- 
lution is written in terms of T and n as 



(33) 
(34) 

(35) 



Note that these solutions form simple generalization 
of the v is) = 1 case of the solutions of Ref. [1] . Also 
note that in the case when p = nT and n is conserved, 
for some choices of the k (T) function our solution 
becomes ill-defined. The criterion that (k (T) T) 
should be positive limits the applicability of solutions 
for the case of conserved n presented here. In the case 
when for some T range -A; (n (T) T) becomes nega- 
tive, the implicit form of Eq. (35) cannot be inverted 
to give a unique T (r) function. Such domains of T in- 
deed might exist in some parametrizations of the lat- 
tice QCD Equation of State around the quark-hadron 
transition temperature (as detailed in the next sec- 
tion, in particularly on Fig. 1). However, even for 
these cases, one can use the solution without con- 
served n presented here as a physically relevant solu- 
tion, since at the transition temperature a conserved 
density n yielding pressure as p = nT cannot be iden- 
tified. 

Let us briefly mention another possibility, when 
k is a function of the pressure p and not that of the 
temperature T. In this case a solution can be written 
up, similarly to the previous ones as 



(36) 
(37) 
(38) 



i.e. almost the same as in Eq. (32), except that here 
the integration variable is the pressure p. If however 
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Figure 1: The temperature dependence of the EoS 
parameter k from Ref. [7] is shown with the solid 
black curve. Note that in the shaded T range (173 
MeV - 230 MeV) (k (T) T) (red dashed line) be- 
comes negative, thus the implicit form of Eq. (35) 
cannot be inverted to give a unique T (r) function. 
Hence we will substitute this k(T) in the hydrody- 
naniic solution shown in Eqs. (30)-(32). 



the pressure can be written as a function of temper- 
ature, i.e. as p(T), an integral-transformation can 
be made with and we get back Eq. (32), so in this 
case these solutions are identical. This solution may 
be used if a n(p) function is given (without relation 
to the temperature) by an arbitrary energy density 
function e(p) = n(p)p. 



5 Application 

Recently a QCD equation of state has been calcu- 
lated by the Budapest- Wuppertal group in Ref. [7]. 
Here (in their Eq. (3.1) and Table 2) they give a 
parametrization of the trace anomaly as a function of 
temperature. Hence the pressure, the energy density 
and finallly the EoS parameter k can be calculated, 
as a function of the temperature. We did this calcu- 
lation, and got the k(T) function as shown in Fig. 1. 
Note however, that in this calculation for some T 
range ^ (k (T) T) becomes negative, as also shown 
in Fig. 1. Hence the implicit form of Eq. (35) can- 
not be inverted to give a unique T (r) function. We 
can still use the solution without conserved number 
density n, presented in Eqs. (30)-(32). 

We utilized the obtained k(T) and calculated the 
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Figure 2: Time dependence of the temperature T(t) 
(normalized with the freeze-out time r and the 
freeze-out temperature To) is shown here. The four 
thin red lines show this dependence in case of con- 
stant k values, while the thicker blue lines show re- 
sults based on the EoS of Ref . [7] . The resulting curve 
slightly depends on the value of To- It is clear how- 
ever, that the temperature fall-off is almost as fast in 
the QCD EoS case as in the case of fixed k = 3, which 
resembles a relativistic ideal gas. This means, that a 
fixed freeze-out temperature (which cannot vary too 
much due to the known quark-hadron transition tem- 
perature) results in a very high initial temperature. 



We have presented the first analytic solutions of the 
equations of relativistic perfect fluid hydrodynamics 
for general temperature dependent speed of sound 
(ie. general Equation of State). They can be seen 
as generalizations of previously known exact solu- 
tions [1]. However, our new solutions are spherically 
symmetric, thus possible generalizations of them def- 
initely are worth exploring: solutions for the non- 
accelerating case and for more general ellipsoidal 
symmetry would be able to analytically explore the 
time evolution of other hadronic observables such as 
the elliptic flow (t> 2 ). 

We have shown how to use our solutions to fully 
utilize the lattice QCD Equation of State for explor- 
ing the initial state of heavy-ion reactions based on 
the reconstructed final state in the Buda-Lund hy- 
drodynamical model. In ^/snn = 200 GeV Au+Au 
collisions, our investigations reveal a very high ini- 
tial temperature consistent with calculations based 
on the measurement spectrum of low momentum di- 
rect photons [8]. If given a temperature-dependent 
direct photon emission function, then this model can 
be used to calculate direct photon spectra to be com- 
pared to measurements, as in Ref. [9], but with a 
realistic Equation of State. 
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time evolution of the temperature of the fireball from 
this solution of relativistic hydrodynamics. The re- 
sult is shown in Fig. 2. Clearly, temperature falls off 
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A The entropy and the temper- 
ature equations 

The fundamental thermodynamical relations con- 
necting e, T, a, p, and any types of rii conserved 
charges and corresponding /i.; chemical potentials are 

e + p = Ta + nj/Zj, (39) 

i 

de = Tda + ^2mdm, (40) 

i 

dp = adT + ^2n i diJ, i . (41) 
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In the case when there are no conserved charges, sim- 
ilar relations hold with all rij and Hi variables omit- 
ted. Substituting these in the (5) energy conserva- 
tion equation, we immediately obtain the continuity 
equation for the entropy density a: 

Tad^ + TvPdpO + ^ m (n^u" + u^m) = 0, 

(42) 

which is, for conserved (or vanishing) ms, equivalent 
to 



d v {ou v ) = 0, 



(43) 



which is the entropy conservation, Eq. (6). 

In the case when there is no conserved n, we can 
substitute the (40) and (41) thermodynamic relations 
for vanishing n in Eq. (5). Using the Eos as e = 
k (T) p, and e + p = Ta and dp = a AT we obtain 
from Eq. (5) the following: 



Ta 



1 c\k 

^ + ^TTdr^ r 



+ na u^d^T = 0, 



(44) 



which is, by using Eq. (14) again, equivalent to Eq. 
(17), as was to be demonstrated. 

Next, we would like to obtain an equation for the 
temperature with our specific Equation of State as 
in Eq. (7) (e = k (T)p), in the case when there is a 
conserved n and p = nT. We can substitute these 
into Eq. (5), and use the (1) continuity equation for 
n to infer that Eq. (5) is equivalent to the following: 



Td^ 



AT 



(kT) u^d^T = 0. 



(45) 



Introducing V by using Eq. (14), we immediately see 
that this is equivalent to Eq. (16) , as stated in the 
text. 

Finally, let us show how the solution for a given 
n(p), described in Eqs. (36)-(38) can be obtained. 
In that case, instead of substituting the temerature 
to Eq. (5), we write up the equation using the n(p) 
function and the relation e = k ■ p, similarly to the 
previous cases: 
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+ 



d v p 



k dn 
p dp ) k + 1 



= 0. 



(46) 



This equation is then solved by the implicit formula 
on the pressure, given in Eq. (38) 
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